Model for hypernucleus production in heavy ion collisions 
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We estimate the production cross sections of hypernuclei in projectile hke fragment (PLF) in 
heavy ion coUisions. The discussed scenario for the formation cross section of A hypernucleus is: (a) 
A particles are produced in the participant region but have a considerable rapidity spread and (b) A 
with rapidity close to that of the PLF and total momentum (in the rest system of PLF) up to Fermi 
motion can then be trapped and produce hypernuclei. The process (a) is considered here within 
Heavy Ion Jet Interacting Generator (HIJING/BB) model and the process (b) in the canonical 
thermodynamic model (CTM). We estimate the production cross sections for a hypernucleus j^'' Z 
where Z — 1,2,2 and 4 for C + C at total nucleon-nucleon center of mass (cm.) energy ^/snn = 3.7 
GeV, and for Ne+Ne and Ar+Ar collisions at ^/snn = 5.0 GeV. By taking into account explicitly 
the impact parameter dependence of the colliding systems, it is found that the cross section is 
different from that predicted by the coalescence model and large discrepancy is obtained for XVLe 
and ^Be hypernuclei. 

PACS numbers: 25.75.-q, 21.80.+a, 25.70Mn, 25.70Pq 



I. INTRODUCTION 

Hypernuclear production in reactions between heavy 
nuclei was first theoretically proposed by Kerman and 
Weiss [l|. They found that relativistic heavy ion col- 
lisions (rHIC) offer the best possibility to create exotic 
finite nuclear system with finite strangeness. However, 
the experimental results have so far been rather limited 
due to short lifetime of hypernuclei |2j-lJ] , which impedes 
their detection. Since the mechanism of the heavy ion in- 
duced hypernuclear reaction is not well understood, this 
field of research attracted mainly theoretical interest, see, 
e.g., Refs. 

The future experimental projects as planned at the 
Facility for Anti-protons and Ion Research (GSI-FAIR, 
Germany) (23 - [27[ and the Nuclotron-based Ion Collider 
Facility (NICA), at Joint Institute for Nuclear Research 
(JINR), Dubna, Russia aim to look for light pro- 
ton/neutron rich exotic hypernuclei, and to extend their 
study to heavier hypernuclei (which can be produced only 
in rHIC), toward protons and neutron drip lines. In addi- 
tion, there are special light hypernuclei of interest, whose 
properties are dictated by their unusual structure (29j . 
e.g., nuclei with the hyperon halo (^H); neutron- rich hy- 
pernuclei; nuclei with an unstable core, where the hy- 
peron is a sort of "glue" ensuring stability (^H, %iie, 
\Re). Recently at RHIC BNL (Relativistic Heavy Ion 
Collider at Brookhaven National Laboratory) the obser- 
vation of (arLti )hy pertritons has been reported in nuclear 
collisions [30] , 31 1 • 



The research program at the SPS CERN-Geneva (Su- 
per Proton Synchroton at European Organization for Nu- 
clear Research) and at RHIC BNL will cover the energy 
range ^/sNN — 5-50 GeV from below up to well above 
the energy of the onset of deconfinement (expected ap- 
proximately at = 6-8 GeV). Moreover, HypHI- 
coUaboration [22, |23| and the project Multi-Purpose De- 
tector (MPD/NICA) [2^ will cover the relevant energies 



ySjvjv = 2.7-9.4 GeV. These experiments aim to measure 
the production of hypernuclei in energetic collisions be- 
tween light nuclei, since their identification via the weak 
decay of the hyperon into pion is much easier for such 
systems. 

In this paper a hybrid model based on participant- 
spectator picture and combined with canonical thermo- 
dynamic (CTM) model is used to determine, in high en- 
ergy heavy-ion collisions, the probability of forming a hy- 
pernucleus and to estimate its production cross section. 
For rHIC in 3-10 GeV energy range the following scenario 
(supported by experiments) applies. For a general impact 
parameter there is a region of violent collision called the 
participating region. In addition, there is a mildly ex- 
cited projectile like fragment (PLF) and also a mildly 
excited target like fragment (TLF). Physics of both PLF 
and TLF are similar for symmetric collisions; here we 
concentrate our analysis on PLF. Because of excitation 
energy (usually characterized by a temperature, T) PLF 
will break up into many fragments and the velocities of 
these fragments are centered around the velocity of the 
projectile. In fixed target experiments they are emit- 
ted in a forward cone and are more easily recognizable. 
In the participating region, apart from original neutrons 
and protons, others particles (pions, kaons, A's, etc.) are 
produced. The produced A particles have an extended 
rapidity range. Those produced in the rapidity range 
close to that of the projectile and having total momenta 
in the PLF frame up to Fermi momentum (ptot < 250 
MeV/c) can be trapped in the PLF and form hypernu- 
clei. The present problem has also been analyzed previ- 
ously (see Refs. 0, H, [ll|, [l2|) at lower ^/SNN energies 
{\/snn — 2.7-3.1 GeV). This work is in a similar spirit 
but at higher cm. energy and uses different models for 
(a) A particles production in the participant region tak- 
ing into consideration impact parameter dependence of 
the specific collision and (b) attachment of the A particle 
to different fragments in the PLF forming what we refer 
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here as composite. It is felt that in view of the future 
experimental activities (22| - [28| the results from alterna- 
tive models, and at different energies, will help to better 
establish the mechanism for formation of hypernuclei. 

Our calculations are performed in two separate stages. 
For a given impact parameter (b) a large number of events 
(10^ - 10^) are generated with HIJING/BB model m or- 
der to obtain the average number of A particles per event 
(< n\{b) >) within appropriate kinematic cuts in ra- 
pidity and momentum. Within the model we can also 
calculate the average number of the nucleons in the PLF, 
< "-PLF > — < 'ni{b) > + < n2{b) >, where ni{b) stands 
for neutrons, and n2{b) for protons. 

We can characterize a produced composite (with and 
without strangeness) by three symbols: a, z and h, where 
a is the atomic mass, z is the charge value of the isotope, 
and h refers to the number of hyperons embedded into the 
nucleus (a, z). For h = we have normal composites. For 
h = 1 we have a hypernucleus with one A. It is also pos- 
sible to include the case h — 2 (i.e., %j<^z). Similarly it is 
possible to include other kinds of hyperons, i.e., S. How- 
ever, this is beyond the scope of the present work. Cal- 
culations using canonical thermodynamic model (CTM) 
[20I 111] are performed in the second stage to estimate 
the average number (< Ua.z.iib) >) of hypernuclei (^z) 
produced at a given temperature T, when the PLF has 
one A, and the average number of nucleons < ripLF > 
=< ni{b) > + < n2{b) >. 

There are two major approximations in our approach: 
(i) the time dependence of source function for hyperons 
and fragments has been neglected; (ii) we neglect also 
secondary hypernuclear processes. The transport calcu- 
lations predict only moderate contributions [21] to the to- 
tal hyper-fragment cross section from indirect coalescence 
through the ttB channel (where B stands for a baryon and 
IT for a pion). With these assumptions, the formation 
cross section of a hypernucleus in the PLF rest system 
can be expressed by: 



aii^Z) = J 2Trbdb < UAib) > < ria^b) > (1) 

The basic outline and main parameters of HIJING/B13 
model are presented in Sec. II. We employ reasonable 
theoretical evaluations for the A particle production with 
specific kinematic cuts. In Sec. Ill we give details of 
CTM calculations. The results and the discussions are 
presented in Sec. IV. The final conclusions are summa- 
rized in Sec. V. 



II. OUTLINE OF HIJING/BB V2.0 MODEL. 

The HIJING ^ and HIJING/BB vl.lO models [slj 
have been used extensively to determine the physical 
properties of the matter produced in rHIC and to study 
particle production. String models describe the collisions 



through the exchange of color or momentum between par- 
tons in the projectile and target. As a consequence of 
these exchanges, these partons become joined by color- 
less objects, which are called string, ropes or flux tubes. 
In HIJING type models 's^] the soft beam jet fragmenta- 
tion is modeled by diquark-quark strings with gluon kinks 
induced by soft gluon radiations. Hard collisions are in- 
cluded within perturbative Quantum Chromodynamics 
(pQCD) computed parton-parton collisions. The mini- 
jets physics is embedded via an eikonal multiple collision 
framework using pQCD to compute the initial and final 
state radiation and hard scattering rates. The cross sec- 
tion for hard parton scatterings is enhanced by a factor 
K=2 in order to simulate high order corrections. The HI- 
JING model was originally designed for hadron-hadron 
interactions. Generalization to hadron-nucleus (p + A) 
and nucleus-nucleus {A + A) collisions is performed |34| 
through the Glauber theory. Besides the Glauber nu- 
clear eikonal extension, shadowing of nuclear parton dis- 
tributions is modeled. In addition, dynamical energy loss 
of the (mini)jets is taken into account through an effec- 
tive dE/dx (approximately 2 Gev/fm per gluon mini-jet). 
These models contain a soft and a hard component, which 
is crucial for their application from FAIR to Large Hadron 
Collider (LHC) energies. 

Unlike conventional di-quark fragmentation imple- 
mented in HIJING model |34l |. a baryon junction allows 
the di-quark to split with the three independent flux lines 
tied together at a junction in HIJING/B13 vl.lO We 
introduce [Sgl a new version (v2.0) of HIJING/BB that 
differs from HIJING/BB vl.lO [s^ in its implementation 
of hypothesized junction anti-junction (JJ) loops. A pic- 
ture of a junction loop is as follows: a color flux line splits 
at some intermediate point into two flux line at one junc- 
tion and then the flux line fuse back into one at a second 
anti-junction somewhere further along the original fiux 
line. Before fragmentation, we compute the probability 
that a junction loop occurs in the string. The probabil- 
ity of such loop is assumed to increase with the number 
of binary interactions that the incident baryon suffers in 
passing through the oncoming nucleus j36| . This number 
depends on the impact parameter and is computed in HI- 
JING using eikonal path through a diffuse nuclear density. 
Moreover, we add an enhanced intrinsic (anti)di-quarkpT 
kick in any standard (q-qq) strings that should contain 
one or multiple JJ loops. This mechanism correspond to 
some degree of collectivity for parton interactions. 

In string fragmentation phenomenology, it has been 
proposed that the strong enhancement of strange particle 
observables require strong color field (SCF) effects [s^l . 
The particle production is large and dominated by pair 
production and the energy density appears to exceed the 
one required for quark gluon plasma (QGP) formation. 
The overall scenario is consistent with particle produc- 
tion from a strong color field (SCF) , formation of a QGP 
and subsequent hadronization. The SCF effects are also 
possible source of novel baryon/hyperon physics. In our 
previous works f38l - l4l| . we explore with HIJING/BB v2.0 
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model calculations, possible dynamical effects associated 
with long-range coherent field (SCF), that may arise in 
rHIC. These analysis point to the need for collective mo- 
tion, together with high initial-state energy densities. 

For a uniform chromoelectric flux tube with field (E) 
the pair production rate ^42] per unit volume for a heavy 
quark is given by: 



r 



47r3 '^^P 



(2) 



where for strange quark Q — s, the current quark mass is 
in the range of TOs = 80 — 190 MeV and the constituent 
quark mass Ms = 350 MeV (for detailed discussion see 
Ref. ;39j). Note that k = \eE\eff ^ \/ C2{A) / C2{F) kq 
is the effective string tension in terms of the vacuum 
string tension kq ~ 1 GeV/fm and C2{A), C-zjF) are 
the second order Casimir operators (see Ref. [23]). A 
measurable rate for spontaneous pair production requires 
strong chromo electric fields" , such that k/tuq > 1 at 
least some of the time. On the average, longitudinal elec- 
tric field "string" models predict for heavier fiavor a very 
suppressed production rate per unit volume 7q via the 
well known Schwinger formula [4^ . since 



gg 
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for Q = s and q = u,d. For a color rope on the other 
hand, if the average string tension value (< k >) in- 
creases from 1.0 GeV/fm to 2.0-3.0 GeV/fm, the rate F 
for strangeness pairs to tunnel through the longitudinal 
field increases (see Refs. (ssl l39j). 

We take into account SCF in our model by an in 
medium effective string tension k > kq, which lead to new 
values for the suppression factors, as well as the new ef- 
fective intrinsic transverse momentum fc^. This includes: 

i) the ratio of production rates of di-quark to quark pairs 
(di-quark suppression factor), 7qq = -P(qqqq)/-P((Zq), 

ii) the ratio of production rates of strange to non- 
strange quark pairs (strangeness suppression factor), 
7s — P{ss)/P{qq), in) the extra suppression associ- 
ated with a diquark containing a strange quark com- 
pared to the normal suppression of strange quark (7s), 
Jus = (-P(usus)/P(udud))/(7s), iv) the suppression of 
spin 1 diquarks relative to spin ones (apart from the fac- 
tor of 3 enhancement of the former based on counting the 
number of spin states), 710, and v) the (anti)quark (cr^' = 

a/k/kq • fq) and (anti) di-quark {cr'^^ — ^/k/kq ■ f ■ o-qq) 
Gaussian width. 

In this paper we extend our study in the framework 
of HIJING/BB v2.0 model ^ to strange particle produc- 
tion at the FAIR and MPD/NICA energy range (-^ snn 
= 2.7-9.4 GeV). The experimental situation has so far 
been rather poor for measurements of strange particles 
at forward rapidities, due to limited acceptances of the 
detectors. The very forward rapidity (closer to projec- 
tile rapidity) are dealt with by models only. Using HI- 
JING/BB v2.0 model we analyze the production of the 



TABLE I: Main parameters used in symmetric ^4-1-^4 collisions. 
The parameters are defined in the text. Set 1 correspond to 
calculations without JJ loops and SCF eff'ects. Set 2 and 3 
include both effects and correspond to different mean values 
of the string tension. 



A + A 


(GeV/fm) 


7qq 


7s 


7us 


710 


ag 
(GeV/c) 


/ 


Set 1 


KO = 1.0 


0.02 


0.30 


0.40 


0.05 


0.350 


1 


Set 2 


K2 = 2.0 


0.14 


0.55 


0.63 


0.12 


0.495 


3 


Set 3 


K3 = 3.0 


0.27 


0.67 


0.74 


0.18 


0.606 


3 



average value per event < n\{b) > with specific kine- 
matic cuts in rapidity (j/pr-oj —0.05) < Sy < (j/proj +0.05), 
(where yproj stands for projectile rapidity) and total mo- 
mentum ptot < 250 MeV/c in the PLF rest system, which 
can leads to the formation of hypernucleus in rHIC. The 
main parameters used in the calculations presented here 
are given in Table HI 

In these calculations we do not consider nuclear effects 
such as shadowing and quenching (3^ , (36j , which are spe- 
cific only for very high energy interactions {^^SNN > 20 
GeV). Within our model we do not include rescattering 
processes which Relativistic Quantum Molecular Dynam- 
ics (RQMD) model calculations 43] show to be negligi- 
bles at the forward rapidity. The physics embedded in 
our model include a "hard"(pQCD) and an underlying 
"soft" part, therefore we can not apply our phenomenol- 
ogy for energies ^/SNN < 3.0 GeV. The detailed analysis 
for (multi)strangeness production from FAIR to SPS and 
RHIC energies within HIJING/B13 v2.0 model will be pre- 
sented elsewhere [43 |. 

For a given impact parameter (b) we calculate with 
HIJING/BB v2.0 model the average number per event < 
n\{b) >, in appropriate kinematic phase space, and the 
average number of nucleons in the PLF , < np]^p(b) >. 
The predictions for transverse momentum (pr) spectra of 
A particles with rapidity cut Sy for symmetric colliding 
systems ^^C + ^^C, ^ONe + ^ONe, and '^^Ai + ^OAr at 
y^SNN — 5 GeV are presented in Fig. [T] for central {b — 
3 fm, left panel) and peripheral (& = 6 fm, right panel). 
The values of the parameters used to obtain the results 
discussed here, are Set 2 from Table HI 

The average value per event < nA(6) > within rapidity 
cut Sy and with total momentum (ptot < 250 MeV/c) 
in the projectile rest system, are given in Table [H] for 
the symmetric colliding system -^^C -I- -^^C, ^"^Ne -I- ^"Ne 
and '*'^Ar + ^''Ar. These numbers are the main input in 
Eq. 1 used to calculate inclusive cross sections for the 
formation of hypernuclei in nucleus-nucleus collisions in 
the PLF rest frame. 

Our estimate for total number of events with one A 
particle produced in the projectile rapidity region in C -I- 
C collisions at yi^vw = 3.7 GeV is 13.05 x 10""*. A com- 
parable value (26.6 x 10"'*) at lower energy (2A GeV in 
laboratory system) , is obtained [l^ with ultrarelativistic 
quantum molecular dynamics (UrQMD) calculations by 
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FIG. 1: (color online). The HIJING/BB v2.0 model predic- 
tions for pt spectra of As particle within rapidity cut Sy at 
y^sjvjv = 5 GeV. For clarity the values for Ne + Ne and C + 
C are multiplied by 0.1 and 0.01 respectively. 



TABLE II: The dependence on the impact parameter (b) of the 
mean number of A's per event (< nA{b) >) obtained within 
kinematic cuts (see Sec. II) in the PLF rest frame for sym- 
metric colliding systems: ^^C + ^^C, ^"Ne + ^°Ne and '*''Ar 
+ *"Ar. 



C + C 
Vsivjv 3.7 GeV 
6 (fm) < n1{b) > xlO* 


Ne -1- Ne 
5.0 GeV 

<<°(6) > xlO* 


Ar -1- Ar 
5.0 GeV 

< n^'(6) > X 


0.0 


4.61 


6.70 


16.4 


1.0 


2.80 


6.00 


16.20 


2.0 


3.01 


5.60 


18.90 


3.0 


1.44 


5.20 


15.80 


4.0 


0.38 


3.40 


11.80 


5.0 


0.42 


2.00 


9.92 


6.0 


0.24 


0.77 


4.26 


7.0 


0.063 


0.42 


2.74 


8.0 


0.088 


0.25 


1.62 


9.0 




0.095 


0.49 


10.0 






0.23 


11.0 






0.13 



assuming a coalescence of the produced A hyperons with 
spectators. However, their rapidity cut is much larger 
iVproj - 0.4) < SyQMD < ivproj + 0.4), than that used 
in our calculations. Therefore, we investigate within HI- 
JING/BB v2.0 model the sensitivity of our predictions to 
different rapidity cuts Sy. An increase of 6y from 0.05 
to Sy = 0.5 results in an increase of a factor 6-8 for the 
average value < nA(fc) >. The sensitivity to the mean 
values of the in medium string tension (k) is less impor- 
tant. Using Set 3 from Table U which correspond to a 
mean value k = 3 GeV/fm, result in a modest increase 
(less that 3-5%) of < nA{b) >, for Ne -I- Ne central 
collisions (6 = 3 fm) at ^/snn — 5 GeV. 

Finally, for the symmetric colliding system -'^^C-I-^^C, 
we investigate also the energy dependence of particle pro- 
duction for central collisions (6 = 3 fm). When ^/snn 
increase from 3.7 GeV to 5 GeV, the number of A parti- 
cles (< n\{b) >) within kinematic cuts discussed above, 
has a modest increase of only 32 %. The model predict a 



more dramatic increase of < n\{b) >, i.e., a factor of 2.75 
when y/SNN increase from 5 GeV to 10 GeV. Since the 
energy range of the onset of deconfinement is expected 
to be at ^/sNN — 6-8 GeV, the above energy range (5- 
10 GeV) should be carefully investigated by the future 
planned experiments. 



III. MODEL FOR HYPERNUCLEUS 
FORMATION IN THE PLF 

The A particles with appropriate rapidity find them- 
self in the PLF's. They thermalises along with the nucle- 
ons. In an event the A hyperon can remain unattached 
to nucleons or can combine with some nucleons to form 
a hypernucleus. There will be also "normal" composites, 
those without strangeness. Assuming equilibration we 
can compute the average numbers of both normal com- 
posites and hypernuclei. Past experience has shown that 
the temperature in the PLF is expected to be in the range 
of 5 to 8 MeV 

The Canonical Thermal Model (CTM) for two kinds of 
particles (neutron and proton) has had long usage [isf . 
The extension to three kinds of particles (neutron, proton 
and A) was discussed by us in Ref. jl^l- In this analysis 
we just give the details necessaries to follow the calcula- 
tions performed here. In CTM we need to calculate the 
partition function Q: 



Q 



Af,Z,H — 



En 



K,../^)""■ 



na,z,h 



(4) 



where Ap is the number of nucleons in the PLF plus one 
(the A particle); Z is the number of protons in the PLF 
and H = I (only one A with appropriate kinematic cuts 
discussed in Sec. II entered the PLF); h = (normal 
composites) or ft, = 1, a hypernucleus; tOa.z.h is the par- 
tition function of one composite with quantum numbers 
a, z and /i; na,z,h is the number of this composite in a 
given channel. The sum in the above equation includes 
an enormous number of channels which satisfy conserva- 
tion of quantum numbers: 



'^hUa^zM = H 



(5) 



It is shown in Ref. [20[ that the average occupation of 
(a, z, h) is given by: 



< na^z,h >- 



1 



{Af,Z,H 



Ja,z,hQAF-a,Z-z,H-h (6) 



The partition functions Q can be calculated using recur- 
sion relations discussed in Ref. f20j . 

The one particle partition function LOa.z.h has two parts: 



{a,z,h) z^nt{a,z,h) 



(7) 
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The kinetic part is given by 



Zk^n{a,z,h)^^{27rMT)^^^ 



(8) 



where V is the freeze-out volume in which thermodynam- 
ics is assumed. For an atomic mass a we take V as 3 times 
the normal nuclear volume { Vq — a/po, with po « 0.15 
fm^''), but then this volume is reduced for Van der Waals 
excluded volume {V = 2Vb). The mass of the composite 
M is: M = {a — h)mn + hm\, where to„ (938 MeV) and 
TOA (1116 MeV) stand for the nucleon mass and A mass 
respectively. 

The internal partition function, Zintio,, z,h) could be 
written as: 



Zintia, z, h) = ^(2si -f l)exp{-f3ei) 



(9) 



where are energy eigenvalues of the composite (a, z, h). 
However, considerable caution is needed in taking the 
sum over the energy states i. Because of Levinson's the- 
orem the sum needs to be modified and cut off (see the 
discussion on Ref. [1^). In the previous work [2^ the 
interest was to investigate gross features and a relative 
production only was discussed. Therefore, a generic for- 
mula for ground state binding energy and excited states 
was used for most nuclei. 

In this analysis we are attempting a more quantitative 
estimate for production cross sections of hyperfragments 
for Z = 1 to 4. The calculated cross sections differ 
significantly if we use the generic formula from Ref. [20] , 
or experimental energies for ground and excited states. 
In the calculations performed here we have taken exper- 
imental values of ground state and excited state energies 
for low mass hypernuclei (a < 11). We include only par- 
ticle stable excited states. The binding energies of ground 



states for composites with strangeness (a, z, 1) are esti- 
mated from tabulated values of and binding energy 
of normal composites (a — 1, z, 0). Data for a < 11 are 
taken from Refs. [1, and references therein. The de- 
tails are: 

^H: only 1 state with spin 1/2. 

^H: ground state with spin and one excited state with 
spin 1 at 1.04 MeV excitation. 

^He: ground state with spin and one excited state 
with spin 1 at 1.15 MeV excitation. 

^He: only ground state with spin 1/2. 

^He: ground state with spin 1. 

^Li: ground state with spin 1. 

^He: ground state with spin 1/2. 

^Li: ground state with spin 1/2, excited state spin 3/2 
at 0.69 MeV, excited state spin 5/2 at 2.05 MeV, excited 
state spin 7/2 at 2.521 MeV, excited state spin 1/2 at 
3.56 MeV. 

^Be: ground state with spin 1/2. 

^Li: ground state with spin 1, one excited state with 
spin 1 at 1.22 MeV excitation 

^Be: ground state with spin 1, one excited state with 
spin 1 at 1.22 MeV excitation. 



^Li: ground state with spin 3/2. 



^Be: ground state with spin 1/2 and excited states 
with spins 5/2 and 3/2 grouped at 2.93 MeV. 

lOg. ground state with spin 1, excited states with spin 
2 at 0.22 MeV, another with spin 2 at 2.47 MeV and a 
spin 3 at 2.70 MeV. 

j^'^Be: experimental binding is used but excitation en- 
ergies and spins are taken to be the same as for \"B by 
appealing to isospin symmetry. 

For heavier hypernuclei (a > 10), a liquid-drop formula 
is used for ground state energy. This formula is taken 
from Ref. All energies are in MeV. 
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-16a -I- a{T)a^'^ + Q.72z'^ / {a}/^) + 25(a ~h- 2zf/(a - h) - IQMh + 21.27/i/(a^^^) 



(10) 
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where a{T) is the surface tension which is a function de- 
pendent on temperature: 



(t{T) = 18 



T2 -I- r2 



5/4 



(11) 



A comparative study of the above binding energy for- 
mula can be found in Ref. ^9|] . This formula also defines 
the drip lines. We include all nuclei within drip lines in 
constructing the partition function. 

In order to calculate Zi„((a, z, h) we use the liquid-drop 
formula for ea,z,h and include also the contribution com- 
ing from the excited states. This results in a multiplica- 
tive factor exp (r(T)Ta/eo), where r(T) = 12/(12+r) is a 
correction term. For a detailed discussion of parameters 



used in these formulae see Ref [48j. As temperature T 
increases, this correction slows down the increase of par- 
tition function Zint{a, z, h) due to excited states. Similar 
correction has been used before [i^ H^, although this 
is not important for the temperature range considered 
here. We note, that our calculations take into consider- 
ation the effects of the long-range Coulomb force in the 
Wigner-Seitz approximation |43|. 

Apart from hypernuclei, we need also to specify the 
partition function w for normal composites (/i = 0). For 
iH, 2h, 3h, 3He, ^He, ^He, ^Li, ^He, ^Li and ^Be we use 
experimental ground state energy and no excited states. 
For atomic mass a > 6 we use the generic formula from 
Eq. 10 with h = and consider also the contributions 
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from the excited states as described above. 

The temperature (T) in the PLF is also an important 
parameter in our calculations and it is expected to have 
no dependence on the beam energy value (in the range of 
3-10 GeV). Many data can be used to estimate the tem- 
perature and the range T = 5 MeV to T = 8 MeV is a rea- 
sonable one in describing these data [13] , . Our results 
show, that predicted cross sections can change consider- 
ably within this range of the input temperatures, since at 
lower temperature heavier hypernuclei are favored at the 
expense of lighter ones. The average occupation num- 
ber < ria,! > = J2z'a^ '^ill f^'"^* drop with a, go through 
a minimum and rise again (see Fig. 1 and Fig. 2 from 
Ref. [13]). This is the well-known U shape. As the tem- 
perature changes, the shape changes very quickly to a 
different one, where the occupation falls monotonically 
with increasing a. The temperature domain of 5 MeV to 
8 MeV is precisely the range where this happens. As a re- 
sults occupations of the heavier constituents in the PLF 
can change dramatically. These suggest that the rela- 
tive populations of lighter hypernuclei could be a good 
"measure" of the temperature in the PLF. Such effects 
are also seen in our results discussed in the Sec. IV (see 
Table IIIip , but the situation there is more complex, be- 
cause of impact parameter dependence and the expected 
fluctuations on the nucleon numbers in the PLF. 

We also note, that the PLF system is finite, and 
in this temperature range grand canonical model re- 
sults differ drastically from those obtained within CTM 
model, which include additional constrains such as par- 
ticles number conservation. The detailed comparison of 
the results obtained using both models can be found in 
Ref. m. 



IV. RESULTS AND DISCUSSIONS 

This work presents the inclusive production cross sec- 
tions of different types of hypernuclei for the symmetric 
coUiding systems ^^C -I- ^^C at y/s^ = 3.7 GeV, ^"Ne 
+ 20Ne and '^^Ai + ^°Ai- at ^i]^ = 5 GeV. We assume 
that A particles with a rapidity cut 6y and total momenta 
less than 250 MeV/c in the rest frame of the projectile 
will thermalises with the nucleons in the PLF and will 
produce hypernuclei. In all three cases studied here the 
following characteristics are common. Integration over 
impact parameter is important. The relevant quantity is 
the product b < nA{b) > in Eq. 1. The peak of this func- 
tion is not too sharp. For small impact parameters {b < A 
fm) the number of nucleons in the PLF is small and thus 
only light hypernuclei can be formed when the hyperon 
A is captured. For large impact parameters (& > 4 fm) 
the number of nucleons in the PLF is greater than in 
the above case, and both light and heavy hypernuclei are 
formed. 

Previous theoretical studies through coalescence model 
d, i, M, m predicted the cross sections of the order of 
few microbarn (p-b). The inclusive cross sections obtained 



TABLE III: Inclusive production cross sections (all values are 
in fib) for different types of hypernuclei for the colliding system 
i2(j _j_ i2(-, = 3.7 GeV. Our results (column 2 and 



3) are compared with p: 
coalescence model [9l4lll] 


revious predictions obtained within 
at ^snn ~ 2.7 GeV. 


Model 
Hypernuclei 


Coalescence CTM 

(T = 5 MeV) (T 


CTM 
= 8 MeV) 


iH 


0, 


,22 


0.89 


3.25 


iH 


0, 


.39 


0.32 


0.71 


iHe 


0, 


.39 


0.34 


0.77 


iHe 


2, 


.58 


3.87 


1.46 


iHe 


0, 


.32 


0.50 


0.17 


iHe 


0, 


.09 


0.0009 


0.004 


iLi 


0, 


.30 


0.56 


0.18 


iLi 


0, 


.24 


26.88 


0.85 


iLi 


0, 


.18 


0.17 


0.13 


iLi 


0, 


.05 


0.00008 


0.0004 


iBe 


0, 


.07 


0.001 


0.005 


iBe 


0, 


.15 


0.18 


0.13 


iBe 


2, 


.48 


22.3 


2.26 


i«Be 


0, 


.33 


0.02 


0.018 



in our model for different types of hypernuclei (as indi- 
cated) for the hght colhding system ^^C -I- ^^C at ^snn 
= 3.7 GeV are given in Table Hill for temperature T = 5 
MeV (second column) and T — 8 MeV (third column) in 
comparison with coalescence model results (at y^s/vw ~ 
2.7 GeV) d, i, [Til, [il]- One knows that the total cross 
section of nucleon-nucleon {N — N) scattering becomes 
nearly constant at aproximately 40 mb for cm. ener- 
gies in the range 2.5-10 GeV, and in this energy region 
only 1 % of the {N — N) scattering produces a baryon 
of strangeness 1. Also we note that there is an energy 
threshold of « 1.6 GeV for A production in an elemen- 
tary process of NN AK N. Therefore, we could con- 
sider this comparison as appropriate, since in the region 
^SNN — 2-4 GeV, the A hyperons production is expected 
to have a modest increase. Moreover, for Ne -|- Ne col- 
lisions, the calculations with coalescence model [^ show 
that the hypernucleus formation cross sections have also 
a modest increase (by 30-50 %), when ^/snn increase 
from 2.04 GeV to 3.1 GeV (see Fig. 6 in Ref. [|). 

Our results are different from those reported within 
the coalescence model [1, [1, [ll[, [l2| which was used often 
to estimate cross sections for hypernucleus formation in 
the PLF. Only for few isotopes (e.g., iHe, iLi, iBe) we 
obtain inclusive cross sections slightly different in com- 
parison with those predicted by the coalescence model. 
It is worthwhile to highlight the important differences. 

To understand the basic idea of the coalescence model, 
let us begin by first considering a PLF without a hy- 
peron. It is excited (usually parametrized by a temper- 
ature) and breaks up (as is well established experimen- 
tally) into many fragments. The distribution of these 

fragments (F) in the PLF can be denoted by [pf)- 
Usually this distribution is not calculated in the coales- 
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cence model for hypernucleus formation although in some 
version of the coalescence this could be attempted. In 
principle, this distribution should be taken from experi- 
ment. In the present problem, we have, in addition, a A 
particle with a momentum distribution. If the velocity 
of the fragment pp/rnF and the velocity of the A parti- 
cle, pt^/mf^ match they can coalesce into a hypernucleus. 
Thus the cross-section is given by the product of two dis- 
tributions and an overlap factor which gives the probabil- 
ity that the A and F becomes the hypernucleus aF (see 
Eq. 15 in Ref. and Eq. 11 in Ref. 0). Note that the 
A plays a very passive role here. Practitioners of the coa- 
lescence model suggested that has to be taken from 
experiment. In the past applications of the coalescence 
model, experimental values of ^3^^ were not used as not 
all of those needed for predicting hypernucleus produc- 
tion were experimentally available. A theoretical model 
parametrisation for ^jt— had to be employed. These pa- 
rameterizations give (wrongly) finite values for inclusive 
cross-sections of ^He and *Be nuclei and hence finite val- 
ues for inclusive cross-sections formation for both ^He 
and A^e. But this is fortuitous, since ^He and ^Be nu- 
clei are not bound systems by themselves [13, [E^l ■ How- 
ever, calculations based on a generalized mass formula for 
non-strange and hypernuclei predict the existence of sev- 
eral bound hypernuclei (e.g., A He and ^Be) whose normal 
cores are unbound (sil . [52j . 
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FIG. 2: (color online). The hypernucleus formation cross sec- 
tions ((JHy) in ^^C + ^^C collisions at ^/snn = 3.7 GeV as 
a function of nuclear fragment mass number Ap- The results 
(open symbols) are shown for H, He, Li and Be isotopes and 
are obtained for T = 5 MeV (left panel) and T = 8 MeV (right 
panel). The lines are only to guide the eyes. 

These remarks lead to some features which are exper- 
imentally verifiable. For example, correct application of 
coalescence principles should gives zero cross section for 
^He and ^Be isotopes, since these require separately a 
^He nucleus and a A or a ^ Be nucleus and a A. The ex- 
tra role the A play as "glue" ensuring stability, can not 
be incorporated in the coalescence model. In contrast in 
CTM model the A plays a more fundamental role than 
in coalescence. It participates in the thermalization. The 
CTM model uses directly the property of aF rather than 

^He and ^I 



expected to have non-zero production cross sections. In 
fact, we predict a large cross section for A^e. The val- 
ues obtained within our phenomenology are also shown 
in Fig. ma) (T = 5 MeV) and in Fig. HJb) (T = 8 MeV). 
Note that the high value reported here for isotope ^Li 
(on dashed lines) could be explained by many low-lying 
bound states of this isotope, which all contribute in the 
calculation of the partition function, Zint- 

In collisions between heavy nuclei we expect an in- 
crease of the production rate of strangeness (see Table 
mi in Sec. II) and of secondary interactions. In this case 
the temporal distribution of the fragment and A particle 
could play an important role and the reduction factors 
of about 0.5 was estimated with coalescence model [l^ 
for 20Ne + ^^Ne at y/s^N = 3.1 GeV. Our results have 
been obtained by neglecting the above time dependence 
for 20Ne -f ^^Ne and ''"Ar + '^''Ar colhsions at ^/sWn = 
5 GeV. Therefore, these values should be considered only 
as upper bounds for formation cross sections of hypernu- 
clei. The predictions obtained within our hybrid model 
are shown through graphs in Fig.[3]and Fig. 2] (left panel, 
T = 5 MeV) and (right panel, T = 8 MeV). 
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FIG. 3: (color online). The hypernucleus formation cross sec- 
tions in ^''Ne + ^''Ne collisions at y^sjviv = 5 GeV. The open 
symbols and the lines have the same meaning as in Fig. [2] 
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that of F and A separately. Thus both ^He and aBc are 



FIG. 4: (color online). The hypernucleus formation cross sec- 
tions in *°Ar -|- *''Ar collisions at y^ijvjv = 5 GeV. The open 
symbols and the lines have the same meaning as in Fig. [5] 

The Ar -t- Ar system is twice as big as Ne -I- Ne (the 
total number of 80 nucleons versus 40). If we consider 
the results obtained for Li isotopes with T — 5 MeV (left 
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panel of Fig. |3] and Fig. S]), the graphs show that the 
hypernuclei formation cross sections are in Ar + Ar sys- 
tem approximately twice than those obtained in Ne -f Ne 
system. In contrast, the production of ^Be is about the 
same in both. One common feature of all three graphs 
(see Fig. [21 Fig. |3l Fig. |4]) is that, the cross section for 
formation of ^Be is large and is an order of magnitude 
higher than those predicted within coalescence model and 
this prediction should be easily verified by future experi- 
ments planed at FAIR and MPD/NICA. 

Finally, we note that the previous experimental data 
for the production cross sections of and hypernu- 
clei 0, [j] have been obtained with very light beams (He, 
Li) at 3.7A GeV (laboratory system) on carbon fixed tar- 
get. We can not study using our approach these light 
colliding systems, since they are out of the limit of ap- 
plicability within our phenomenology. To compare our 
predictions for formation cross sections of hypernuclei, 
data obtained with heavier beams (targets) and at higher 
energies are needed. 

V. CONCLUSIONS 

We have estimated the inclusive production cross- 
sections of selected light hypernuclei in the symmetric 
colliding systems (C -I- C, Ne -I- Ne, Ar -I- Ar) using a 
hybrid type model. The calculations are performed con- 
sidering a product of two factors in impact parameter 
space (see Eq. 1). The first factor the average value per 
event, < nf^{h) >, estimates the production of A particles 
in the participant region within a limiting phase space. 
Imposing reasonable kinematic cuts these particles could 



be captured in PLF. The second factor, < na^z,i{b) >, 
describes how these hyperons distribute them-self among 
the different hypernuclei within CTM approach. The 
absolute values of cross-sections are basically dependent 
upon the first part; the relative values are fixed by the 
canonical thermodynamic model (CTM). When data will 
become available, these can separately test the validity 
and accuracy of each part. 

In comparison with previous predictions within coales- 
cence models, our results show a striking difference espe- 
cially for ^He and ^Be hypernuclei. A correct application 
of the coalescence principles should gives zero production 
cross sections for these isotopes. In contrast, our hy- 
brid model predict a formation cross sections which are 
significantly larger than those reported with coalescence 
calculations. 

Moreover, our approach for the production of hyper- 
nuclei in rHIC also applies with appropriate changes in 
rapidity and momentum cuts, to asymmetric colliding 
systems. However, to better test our predictions, exper- 
imental data obtained with heavier beams (target) and 
in the energy region of interest ^^snn = 5-10 GeV are 
required. 
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